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Abstract. We study Maxwell's equations on a 4-manifold where the elec- 
tromagnetic medium is described by an antisymmetric (^j-tensor k. In this 
setting, the Tamm-Rubilar tensor density determines a polynomial surface of 
fourth order in each cotangent space. This surfa<;e is called the Presnel surfa<;e 
and acts as a generalisation of the light-cone determined by a Lorentz metric; 
the Presnel surfa<;e parameterises electromagnetic wave-speed as a function of 
direction. Favaro and Bergamin have recently proven that if k has only a 
principal part and if the Fresnel surface of K coincides with the light cone for 
a Lorentz metric g, then k is proportional to the Hodge star operator of g. 
That is, under additional assumptions, the Fresnel surface of k determines the 
conformal class of k. The purpose of this paper is twofold. First, we provide 
a new proof of this result using Grobner bases. Second, we describe a number 
of cases where the Fresnel surface does not determine the conformal class of 
the original (j) -tensor k. For example, if k is invertible we show that k ajid 
have the same Fresnel surfaces. 



1. Introduction 

The purpose of this work is to study properties of propagating electromagnetic fields 
in linear medium. We will work in a relativistic setting where Maxwell's equations 
arc written on a 4-manifold and the electromagnetic medium is represented by an 
antisymmetric (2) -tensor k. Pointwise, such medium is determined by 36 param- 
eters. To understand the propagation of an electromagnetic wave in this setting, 
the key object is the Fresnel surface, which can be seen a generalisation of the 
light-cone [Rub02, HO03, PSW09]. For a Lorentz metric, the light-cone is always a 
polynomial surface of second order in each cotangent space. The Presnel surface, in 
turn, is a polynomial surface of fourth order. For example, the Fresnel surface can 
be the union of two light-cones. This allows the Fresnel surface to model propaga- 
tion also in birefringent medium. That is, in medium where differently polarised 
electromagnetic waves can propagate with different wave speeds. 

The Fresnel surface is determined by the Tamm-Rubilar tensor density which is 
a symmetric (g) -tensor density, which, in turn, is determined by the medium (2)- 
tensor k. This dependence is illustrated in the diagram below: 

Medium k — )• Tamm-Rubilar tensor density Fresnel surface. 

In Lorentz geometry, wc know the the light cone of a Lorentz metric g uniquely de- 
termine 5 up to a conformal factor [Ehr91]. In this work we will study the analogue 
relation between a general electromagnetic medium tensor k and its Presnel sur- 
face; Can one reconstruct an electromagnetic medium from its Fresnel surface? In 
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general, a unique reconstruction is not possible. For example, the Fresnel surface is 
invariant under a conformal change in the medium. Hence the Fresnel surface can, 
at best, determine k up to a conformal factor. One would then like to understand 

the following question: 

Question 1.1. Under what assumptions does the Fresnel surface at a point p G N 
determine the electromagnetic medium k\p up to a conformal factor? 

In terms of physics, Question 1.1 asks when we can reconstruct k\p (up to a con- 
formal factor) using only wavespeed information about the medium at p. A proper 
understanding of this question is not only of theoretical interest, but also of interest 
in engineering applications like electromagnetic tomography. Question 1.1 is also 
similar is spirit to a question in general relativity, where one would like to under- 
stand when the the conformal class of a Lorentz metric can be determined by the 
five dimensional manifold of nuU-geodesics [Low05]. 

Favaro and Bergamin have recently proven the following result of positive nature 
[FBll]: If K has only a principal part and if the Fresnel surface of k coincides with 
the light cone for a Lorentz metric g, then k is proportional to the Hodge star 
operator of g. That is, in a restricted class of medium, the Fresnel surface of k 
determines the conformal class of k. An important corollary is the following: If k 
has only a principal part and its Frcisnel surface coincides with the light cone for 
a Lorentz metric, then k satisfies the closure condition = —/Id for a function 
f:N^ (0, oo). This resolves a conjecture on whether the closure condition char- 
acterises non-birefringent medium in skewon-free medium [OFROO, HO03]. That 
the closure condition is sufficient, was already proven in [OH99, OFROO], but before 
[FBll] sufficiency was only known under additional assumptions; a proof assum- 
ing that ~ (sec Section 2.4 for the definition of in terms on k) is given in 
[OFROO], and a proof in a special class of non-linear medium is given in [OR02]. 
For additional positive results to Question 1.1, see [LH04, Iti05, SWWIO, FBll]. 

The main contribution of this paper is twofold. First, we give a new proof of the 
result quoted above from [FBll]. This is formulated as implication (Hi) =^ (ii) 
in Theorem 4.1. While the original proof in [FBll] relies on the classification of 
skewon-free (2) -tensors into 23 normal forms by SchuUer, Witte, and Wohlfarth 
[SWWIO], we will use Grobner bases to prove Theorem 4.1. Essentially, Grobner 
bases is a computer algebra technique for simplifying a system polynomial equations 
without changing the solution set. See Appendix A. 

The second contribution of this paper is given in Section 5 which contains a number 

of cases, where the Fresnel surface does not determine k. In Theorem 5.1 (iv) we 
show that if k is invertible, then k and have the same Fresnel surfaces. Also, in 
Example 5.3 we construct a k with complex coefficients on M^. At each p e M*, this 
medium is determined by one arbitrary complex number, and hence the medium 
can depend on both time and space. However, at each point, the Fresnel surface of k 
coincides with the usual light cone of the flat Minkowski metric g = diag(— 1, 1, 1, 1). 

The paper is organised as follows. In Section 2 we review Maxwell's equations and 
linear electromagnetic medium on a 4-manifold. In Section 3 we describe how the 
Tamm-Rubilar tensor density and Fresnel surface is related to wave propagation. 
To derive these objects we use the approach of geometric optics. As described 
in Section 3, this can be seen as a step towards a rclativistic theory of electro- 
magnetic Gaussian beams (if such a theory exists). In general, Gaussian beams 
is an asymptotic technique for studying propagation of waves in hyperbolic sys- 
tems. These solutions behave as wave packets; at each time instant, the entire 
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energy of the solution is concentrated around one point in space. When time moves 
forward, the beam propagates along a curve, but always retains its shape of a Gauss- 
ian bell curve. Electromagnetic Gaussian beams are also known as quasi-photons 
[Kac02, Kac04, Kac05, Dah06]. For the wave equation, see [Ral82, KKLOl]. For 
the history of Gaussian beams, see [Ral82, Pop02]. In Section 4 we prove the main 
result Theorem 4.1, and in Section 5 we describe a number of cases where Question 
1.1 has a negative answer. 

This paper relies on a number of computations done with computer algebra. Further 
information about these can be found on the author's homepage. 



2. Maxwell's equations 



By a manifold M we mean a second countable topological Hausdorff space that 
is locally homeomorphic to with C°°-smooth transition maps. All objects are 
assumed to be smooth where defined. Let TM and T*M be the tangent and 
cotangent bundles, respectively, and for fc > 1, let A'^(M) be the set of p-covectors, 
so that A^{N) = T*N. Let ilf (Af) be (j)-tcnsors that are antisymmetric in their 
k upper indices and I lower indices. In particular, let U'^{M) be the set of Re- 
forms. Let also X{M) be the set of vector fields, and let C°°(M) be the set of 
fimctions. By f7'^(M) x K we denote the set of fc-forms that depend smoothly on 
a parameter teR. By T{M,C), T*{M,C), Ap{M,C), nf{M,C) and X(M, C) we 
denote the complexification of the above spaces where component may also take 
complex values. Smooth complex valued fimctions are denoted by C°°(M, C). The 
Einstein summing convention is used throughout. When writing tensors in local 
coordinates we assume that the components satisfy the same symmetries as the 
tensor. 

We will use differential forms to write Maxwell's equations. On a 3-manifold M, 
Maxwell equations then read [BH96, HO03] 

(1) = -f , 

(2) dH = f + J, 

(3) dD = p, 

(4) dB = 0, 

for field quantities E,H (M) x K, D, S e (-^,^) ^ M and sources J e fi^ ^ 
M and p G f2^(M) x M. Let us emphasise that equations (l)-(4) are completely 
differential-topological and do not depend on any additional structure. (To be 
precise, the exterior derivative does depend on the smooth structure of M . However, 
for a manifold M of dimension 1, 2, 3 one can show that all smooth structures for 
M are diffeomorphic. For higher dimensions the analogue result is not true. Even 
for there are uncountably many non-diffeomorphic smooth structures [Sco05, 
p. 255].) 



2.1. Maxwell's equations on a 4-manifold. Suppose E,D,B,H are time de- 
pendent forms E,H & Vl^{M) x M and £), B e ^^'{M) x M and iV is the 4-manifold 
N = WxM. Then we can define forms F,G & VL^{N) and j e ^^{N), 

(5) F = B + EAdt, 

(6) G = D-HAdt, 

(7) j = p- J A dt. 
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Now fields E, D, B, H solve Maxwell's equations equations (l)-(4) if and only if 

(8) dF = 0, 

(9) dG = j, 

where d is the exterior derivative on N. More generally, if A'^ is a 4-manifold and 
F,G,j are forms F,G € CP{N) and j e Q^{N) wc say that F,G solve Maxwell's 
equations (for a source j) when equations (8)-(9) hold. By an electromagnetic 
medium on N we mean a map 

We then say that 2-forms F,G ^ fl^ (N) solve Maxwell 's equations in medium k if 
F and G satisfy equations (8)-(9) and 

(10) G = k{F). 

Equation (10) is known as the constitutive equation. If k is invertible, it follows 

that one can eliminate half of the free variables in Maxwell's equations (8)-(9). We 
assume that k is linear and local so that we can represent k by an antisymmetric 
(2) -tensor k € Q2{N). If in coordinates {a;'}?^Q for A'' we have 

li4i,dx^ O dx™ O ^ ® ^ 
2 9a;* dx^ 



(11) K = -ni^dx' <g) dx" 



and F = Fijdx^ dx^ and G = Gijdx^ dx^ , then constitutive equation (10) reads 

1 

2' 



(12) Gij = -K^^Fr. 



2.2. Decomposition of electromagnetic medium. Let iV be a 4-manifold. 
Then at each point on N, a general antisymmetric (2) -tensor depends on 36 pa- 
rameters. Such tensors canonically decompose into three linear subspaces. The 
motivation for this decomposition is that different components in the decomposi- 
tion enter in different parts of electromagnetics. See [HO03, Section D.1.3]. The 
below formulation is taken from [Dah09]. 

If K e fl2{N) we define the trace of k as the smooth function N given by 

trace k = - n^- 
2 

when K is locally given by equation (11). Writing Id as in equation (11) gives 
Idti = 5i6i - 5i5l, so trace Id = 6 when dim A/" = 4. 

Proposition 2.1 (Decomposition of a (2)-tensors). Let N be a A-mMnifold, and let 

Z = {KGfll{N):uAK{v) = K{u)Avforallu,vGn'^{N), 
tracer = 0}, 

W = {k gQ,1{N) : uAk{v) = -k{u) Av for allu,v eSf{N)} 
= {k e 0.1{N) : uAk{v) = -k{u) Av for allu,v e ^'^{N), 
trace k = 0}, 
U = {flA&nl{N):f€C^{N)}. 

Then 

(13) Oi(iV) = Z ® W ® U, 

and pointwise, dimZ = 20, dimW^ = 15 and dimll = 1. 
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If we write a k S ^2{^) as 



(1)^ + (2)^ + (3)^ 



with (^^K e Z, ^^^K e W, ^^^K e U, then we say that ^^^k is the principal part, ^^^k is 
the skewon part, ^^^k is the axion part of k. 

2.3. The Hodge star operator. By a pseudo-Riemann metric on a manifold M 
we mean a symmetric real (2) -tensor g that is non-degenerate. If M is not connected 
we also assume that g has constant signature. If g is positive definite, we say that 
5 is a Riemann metric. By jj and b we denote the isomorphisms tt: T*M — )• TM 

and b: TM — >■ T*M. By M-linearity we extend g, ft and b to complex arguments. 
Moreover, we extend g also to covectors by setting g{^,ri) = g{^^,ri^) when € 
AliN,C). 

Suppose g is a pseudo-Riemann metric on a orientable manifold M with n = 
dim M > 2. For p E {0, . . . , n}, the Hodge star operator * is the map * : CF{M) 
n"-p{M) defined as [AMR88, p. 413] 

*{dx^' A • • • A dx^") = ^^^^5'^'^ • • • g'''^''ei,...i^i^^,...ijx'''+' A-'-Adx^", 

[n — p)\ ^ 

where x* are local coordinates in an oriented atlas, g = g^dx^^dx^ , det g = det gij, 
g^^ is the ijth entry of (gij)~^, and ei-^-.-i^ is the Levi-Civita permAitation symbol. 
We treat £;i...!„ as a purely combinatorial object (and not as a tensor density). We 
also define e'l'"'" = ei^...i^. 

If 5 is a pseudo-Riemann metric on an oriented 4-manifold N, then the Hodge star 
operator for g induces a (2) -tensor k = *g g Vl\{N). If k is written as in equation 
(11) for local coordinates x"^ then 

(14) 4i = V\9\r9''ea,rs. 

Proposition 2.2. Suppose g is a pseudo-Riemann metric on an orientable 4- 
manifold N. Then *g defines a (^-tensor with only a principal part. 

Proof. Let k be the (2)-tensor induced by *g. Then u A k{v) = k{u) A v for all 
u,v € n'^{N) [AMR88, p. 412]. By Theorem 2.1 it therefore suffices to prove that 
trace K = 0. Let us fix p e TV and let x^ are local coordinates for A'' such that 
g\p is diagonal. If k is written as in equation (11) then equation (14) implies that 
tracer = i/v^j = since is diagonal and Sij^i is non-zero only when ijkl are 
distinct. □ 



A pseudo-Riemann metric g is a Lorentz metric if M is 4-dimensional and g has 

signature (H ) or ( — h -I— 1-). For a Lorentz metric, we define the null cone 

at p as the set G Ap(M, M) : g{S,,£,) = 0}. Usually, the null cone is defined as 
a subset in the tangent bundle. The motivation for treating the null-cone in the 
cotangent bundle is given by equation (37). 

The next theorem shows that the conformal class of a Lorentz metric g can be 
represented either using the (2) -tensor *g or the null cone of g. 

Theorem 2.3. Suppose g, h are Lorentz metrics on an orientable A-manifold N . 

Then the following are equivalent: 

(i) There exists a non-vanishing function A G C^{N) such that h = Xg. 

(ii) *g = *h, where *g and *h ctre the (2) -Censors defined by g and h, respec- 
tively. 



6 



DAHL 



(Hi) g and h have the same null cones. 

Proof. Implications (i) (ii) and (i) (Hi) are clear. Implication (ii) ^ (i) 
proven in [DKS89, Theorem 1], and implication (Hi) (i) is proven in [Elir91, 
Theorem 3]. See also [MS08]. □ 

2.4. Decomposition of k into four 3x3 matrices. Suppose (.x'^, x^, .t^, x^) are 
local coordinates for = M x M such that is the coordinate for M and {x^ ,x^,x'^) 
are coordinates for M. If forms F, G are given by equations (5)-(6), then 

for all i,j = 1, 2, 3 and equation (12) then reads 

(15) Hi 
(16) 

where i,j = 1, 2, 3 and r, s are summed over 1, 2, 3. 

Next we show that in coordinates (.t", x^, x'^,x^) the tensor k is represented by four 
3 X 3-matriccs. To do this, let * is the Hodge star operator induced by the Euclidean 
metric on x^,x'^,x^ so that *dx' = i Eq,6=i e'^^dx'' A dxK Thus B = *dx' 
where = ^e^^'^Bj^ and B^n = ^imnB^ ■ In the same way we define D^,D^,D^. 
Now components D' and B^ represent 2-forms D and B in the basis {*rfa;*}f^i, and 
by equations (15)-(16), 

(17) Hi = ^\{-Er) + ^riB\ 

(18) D' = J^"{-Er) + ^r'B\ 

where i € {1, 2, 3}, r is summed over 1, 2, 3, and 

c^r . _ rO ^ . _ _1 ab ^ri _ _\jab rO « _ Ip. tab mn 

Here r index rows and i index columns in 3 x 3 matrices s^, £§, Ql. Inverting the 
relations gives 

"■Oi ~ ® J> "'Or ~ ^krt "'rs ~ <^krs^ ; "'rs "~ '^krs^ -^l > 

where i, r, s e {1,2, 3} and fc, Z are summed over 1, 2, 3. 

The above matrices jz/, "^^j coincide with the matrices , ^ , defined in 
[HO03, Section D.1.6] and [Rub02]. Since these matrices are only part of tensor 
they do not transform in a simple way under a general coordinate transformation 
in AT (see equations D.5.28-D.5.30 in [HO03]). However, if {a;*}f=o and {P}f=o 
overlapping coordinates such that 



^'(a;^a;^a;^), ie {1,2,3}. 
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Then we have transformation rules 
(19) = 

o:t} dx°- 



(20) SSri = det — 



(21) = det 

(22) %^ = 



9a;" / ' ° dx^ dx^ ' 
dx'^ \ dx'^ dx^ 



dx^ J dx'^ dx^ ' 

9P 9a;° 



° dx^ dx"^ ' 

If '^^^K = then Proposition 2.1 imphcs that k is pointwisc determined by 21 coeffi- 
cients. The next proposition shows that these coefficients can pointwise be reduced 
to 18 when the coordinates are chosen suitably. 

Proposition 2.4. Suppose N is a A-manifold and k G i^i^)- Then 

(i) K has no skewon component if and only if locally 

where is the matrix transpose, and s^, ^ are defined as above. 

(ii) Letp G N. If k has no skewon component, then there are local coordinates 
around p such that is diagonal at p. 

(Hi) Let p G N. If K has no skewon component and g is a Lorentz metric on 
N there are local coordinates around p such that is diagonal and for 
some k e {±1} we have g\p = A;diag(— 1, 1, 1, 1). 



Proof. Part (i) follows by [HO03, Equation D. 1.100]. Since we can always introduce 

a Lorentz metric in local coordinates for N, part (ii) will follow from part (Hi). 
For part (Hi), let be coordinates around p such that g\p = fc diag(— 1, 1, 1, 1) 
for k e {±1}- By (i), matrix S2/\p is symmetric, so we can find an orthogonal 
3x3 matrix P = (P* such that Ps^P^ is diagonal and detP = 1. A suitable 
coordinate system is given by = x° and P = Y^j=i j^^ ■ ^ 



3. Geometric optics solutions 



Let K G f2|(A^) on a 4-manifold N , and let F and G be asymptotic sums 

where P > is a constant, $ G C°°(A^,C) and Ak,Bk € n'^{N,C). Substituting 
F and G into the sourceless Maxwell equations and differentiating termwise shows 
that F and G form an asymptotic solution provided that 



(24) d^AAo = 0, 

(25) d^ABo = 0, 

(26) Bk = nAk, 

(27) d(^> A Ak+i + dAk = 0, 

(28) d^ABk+i+dBk = 0, = 0, 1, . . . . 



In equation (26) we treat /t as a linear map k: Cl^{N,C) Q^{N,C). In equation 
(23) function $ is called a phase function, and forms Ak,Bk are called amplitudes. 
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We will assume that Ini<l> > 0, so that F and G remain bounded even if we take 
P-S- DO. 

Lemma 3.1. Suppose N is a smooth manifold, and let q be a 1-form q G fl^{N,C) 
that is nowhere zero. 

(i) If q /\ A = for some A G Q''{N,C) where k > 1, then there exists a 
{k - l)-form a € n''-'^{N,C) such that A = qAa. 

(ii) If q A a = q A a' for some a, a' G Q,^{N,C), then a = a' + fq for some 
feC°°{N,C). 

Proof. Let tt be the isomorphism T*N ^ TN induced by an auxiliary (positive 
definite) Ricmann metric on N, and let || • || be the induced norms on TN and 
T*N. Let also q = a + i(3, where a = Req and (3 = Imq. Then vector field 
X e XiN, C) given by 

_ a' - iP^ 

NF^PF 

satisfies q{X) = 1. Contracting qAA = OhyX gives part (i). Part (ii) follows by 
taking A = a — a' in part (i). □ 

In this work we will only analyse the leading amplitudes Aq and Bq. However, 
since Bq = k{Ao), it suffices to study Aq in more detail. Let us assume that Aq 
and Bq solve equations (24)-(26). Then Lemma 3.1 (i) implies that there exists a 
1-form oo e Sl^{N,C) such that Aq = A ao, whence 

(29) d^AK{d^Aao) = 0. 

For N = M. X M where M is a 3-manifold and for special choices for k, $ and 

amplitudes A^. , B^. , equation (23) define an electromagnetic Gaussian beam (see 
Section 1). In this setting, $|p and d$|p are both real when p is at a centre of 
a Gaussian beam. With the above as motivation we will hereafter only study 
equation (29) at a point p £ N where is real. From equation (23) we then see 
that d^\p is the direction of most rapid oscillation (or direction of propagation) 
for F. Since Aq = d^ A ao, the 1-form oo. in turn, determines the polarisation of 
the solution in equation (23). Equation (29) is thus a condition that constrains 
possible polarisations once the direction of propagation is known. Since equation 

(29) is a linear in ao, we may study the dimension of the the solution space for 
oq. To do this, let ^ e ^p{N) for some p G N and for ^ let be the linear map 
L^: Ai(iV)^A3(iV), 

(30) L^{a) = CAK(^Aa), aeAl{N). 

We have ^ e kerL^. For all ^ S ApA^\{0} we can then find a (non- unique) vector 
subspace C A^iV such that 

(31) keri| = ® span^. 

Let £^ — d^\p be nonzero. Then Vj\{0} parameterises possible ao that solve equation 
(29) and for which Aq = d^ A Oq is nonzero. For a general k G f^iC-^) ^ ^ 
A^(iV)\{0} wc can have dim E {0, 1, 2, 3}: Proposition 3.7 will show that dim 
can be or 2, Example 3.8 shows that dimV^ can be 1, and the next proposition 
characterise k\p when dimVj = 3 for all ^ e Ap(A'')\{0}. 

Proposition 3.2. Let k G i^N) on a A-manifold N and let p G N. Then the 
following are equivalent: 
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(i) k\p is of axion type. 

(ii) dimVj = 3 for all (, e Ai(7V)\{0}. 

Proof. Implication (i) => (ii) is clear. For the converse direction suppose that (ii) 
holds and a;' are local coordinates around p. It follows that 

CAeAK,(eAa) = 0, a,C,CG Ai(iV). 

If locally ^ = £,idx'^\p then £,i£,jK^ab^^ — ^- Differentiating with respect to and 
gives 

^cr -.dsab , ^dr -csab ri 

'^ab^ "T l^ab^ — ^■ 

With computer algebra it follows that k = | tracer Id and (i) follows. □ 

3.1. Fresnel surface. Let k G f^lC-^) on a 4-manifold N. If k is locally given by 
equation (11) in coordinates a;', let 

C^ijkl _ J^^aiaa aai a^j bib^bak ^b^bibal - 

~ 48 *'l''2 b3&4 6566 ^01020304- 

In overlapping coordinates {S*}, these coefficients transform as 

,^abcd dx' dx^ ds* dx'' 
dx^ ) ° dx'' dx'> dx" dx'^ ' 

Thus components ^^q'''^' define a tensor density on N of weight 1. The Tam,m- 
Rubilar tensor density [Rub02, HO03] is the symmetric part of % and we denote 
this tensor density by 'S^. In coordinates, ^^'-J'^' = where parenthesis indicate 

that indices ijkl arc symmetrised with scaling 1/4!. Using tensor density the 
Fresnel surface at a point p Cz N is defined as 

(33) Fp = U e Ai(7V) : ^^^■'"'='^,^,^6 = 0}. 

By equation (32), the definition of Fp does not depend on local coordinates. Let 
F be the disjoint union of all Fresnel surfaces, F = UpeAf ^p- indicate that Fp 
and F depend on k we also write Fp{K) and F{k). 

If ^ e Fp then e Fp for all A G M. In particular € Fp for each p e N. When 
5f|p is non-zero, equation (33) shows that Fp is a fourth order surface in Ap(A'^), so 
Fp may contain non-smooth self intersections. 

Theorem 3.3. Suppose N is a 4-manifold and k G ^\{N). If S. G ^p(-^) 
non-zero, then the following are equivalent: 

(i) dimVj > 1 where are defined as in equation (31). 

(ii) ^ belongs to the Fresnel surface Fp c K^{N). 

Proof. Let {a;*}|^o be coordinates around p such that dx^\p = ^ and let * be the 
Hodge star operator induced by the Euclidean Riemann metric Qij = 5ij in these 
coordinates. Let P: A^(7V) A^(A'') be the map F = 2 * oL^. Then locally 

3 

F(a) = Y^aiS^^'^^^^dx^, 

where a = aidx'^\p and k^;, are defined as in equation (11). It follows that in the 
basis {rfa;'|p}f^O' the map P is represented by the 4x4 matrix diag(0, Q), where 



(32) 



det 



10 



DAHL 



Q is the 3 X 3 matrix Q"-' = eO'^^-^kO*, i,j e {1, 2, 3}. Now dim > 1 is equivalent 
with dim ker P >2 which is equivalent with det (5 = 0. Writing out det (3 = using 

detQ = yabce^JkQ'''Q'"Q''' 

gives '^^^-'''''^iij^k^.i = 0. We omit the proof of the last step which can be found in 
[Rub02] and [HO03, p. 267 - 268]. □ 

Suppose 5 is a pseudo-Riemann metric on an orientable 4-manifold A'^ and k G 
f2|(-/V). Then g and k define a symmetric (jj) -tensor on A'' by 

(34) = ^ ^-^-^'-^ _ a _ a _ a 



^1 det g\ dx^ dx^ dx'' dx^ 



where ^^'-J'^' are local components of the Tamm-Rubilar tensor density for k, and 
a;' are coordinates in an oriented atlas for A''. 



A key property of symmetric (q) -tensors is that they arc completely determined by 
their values on the diagonal [Muj06, PSW09]. For symmetric (Q)-tensors on a 4- 
manifold (like J^g,^), the precise statement is contained in the following polarisation 
identity. 

Proposition 3.4. Suppose L is a symmetric {^-tensor on a A-manifold N. If 
xi,X2,xs,Xi e Ap(A'') then 

1 ^ ^ 

L{xi,X2,X3,X4,) = ^ ^ 91620364 L(^eiXi,...,^eiXi). 

eie{±i} i=i i=i 

3.2. Electromagnetic medium induced by a Hodge star operator. In Propo- 
sition 2.2 wc saw that a pscudo-Ricmann metric on a 4-manifold induces a (2) -tensor 
K with only a principal part. The next example shows how standard isotropic elec- 
tromagnetic medium can be modelled using a Lorentz metric on R^. 

Example 3.5. On A' = R x ffi^ let «; be the (2) -tensor determined 3x3 matrices 

S2f = -eld, ^ = ir^ld, "^ = ^ = 0, 

where e, /x: — (0, 00). Then constitutive equations (17)-(18) are equivalent with 
the isotropic constitutive equations 

(35) D = e*oE, 

(36) B = n*oH, 

where e is the permittivity and fi is the permeability of the medium and *o is the 
Hodge star operator induced by the Euclidean metric on M.^. If k is the (2) -tensor 

defined as k = ^J~^*g where g is the Lorentz metric g = diag(— 1, 1, 1), then 

equations (35)-(36) are equivalent with equation (10). □ 

The next proposition shows that if 5 is a pseudo-Riemann metric with signature 

(+ + -|— 1-) or ( ) then the medium with k = *g has no asymptotic solutions. 

That is, if d'^\p is non-zero, then equation (29) implies that Ao|p = 0. The propo- 
sition also shows that if k = *g for an indefinite metric g, then Aq can be non-zero 
only when is a null caveator, that is, when g{d^\p,d^\p) = 0. 

Let sgn: K {—1, +1} be the sign function, sgnx = —1 for a; < 0, sgna; = for 
a; = and sgna; = 1 for a; > 0. 



DETERMINATION OF ELECTROMAGNETIC MEDIUM FROM THE FRESNEL SURFACE 11 



Proposition 3.6. Let g and h be pseudo-Riemann metrics on N on an orientable 
A-manifold N. Then 



Proof. Let ^^^J*^' be components for the Tamm-Rubilar tensor density for *g. Com- 
puter algebra then gives 



We know that a general plane wave in homogeneous isotropic medium in can 
be written as a sum of two circularly polarised plane waves with opposite handed- 
ness. The Bohren decomposition generalise this classical result to electromagnetic 
fields in homogeneous isotropic chiral medium [LSTV94]. The Moses decomposi- 
tion, or helicity decomposition, further generalise this decomposition to arbitrary 
vector fields on M^. For a decomposition of Maxwell's equations using this last 
decomposition, see [Mos71, Dah04]. In all of these cases, an electromagnetic wave 
can be polarised in two different ways. Part (i) in the next proposition shows that 
this is also the case for asymptotic solutions as defined above when the medium is 
given by the Hodge star operator of a indefinite metric. 

Proposition 3.7. Let N he an orientable A- dimensional manifold, and let k € 
nl{N) the (2) -tens or K = *g induced by a pseudo-Riemann metric g on N . 

(i) If ^ & A^{N) is non-zero, and is as in equation (31), then 



where = {a & A^{N) : g{a,S^) = 0}. Thus, for any choice of in 



Proof. Let p be the basepoint of ^ and let {x'^}f^Q are local coordinates for N 
around p such that g = gijdx^ ® dx^ and gij \p is diagonal with entries ±1. We know 
that = *l = {-ly Id, where a is the index of g [AMR88, p. 412]. If a e A^(Ar), 
equations (30) and (14) imply that 





^^''"UbUd = sgn(det5)N^d^(ff(C,0)' 
where ^ = $,adx°' and the claim follows by equation (34). 



□ 




equation (31) we have C 




For part (i), equations (37) and (31) imply that dimV^ = dimkerif — 1 where 
H is the 4x4 matrix with entries H'^^ . Let cr{H) denote the spectrum of H 
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with eigenvalues repeated according to their algebraic multiplicity. With computer 
algebra we find that 

where Ci € {±1} are constants that depend only the signature of g. Now part 
(i) follows by Proposition 3.6 and since algebraic and geometric multiplicity of an 
eigenvalue coincide for symmetric matrices [Sza02, p. 260]. For part (ii), equality 
kerLj = follows from the local representation of in equation (37). □ 

The next example shows that the case dimVj = 1 is possible in equation (31). The 
medium defined by equations (38) is called a biaxial crystal [BE80, Section 15.3.3]. 

Example 3.8. On Af = M x M?, let k G Vll{M) be defined by 

(38) = -diag(l,2,3), ^ = Id, '^ = ^ = 0. 
Then the Fresnel equation reads 

(39) 6eo'-eo(5ei' + 8e2' + 9^l) + (e? + el + e;3)(e? + 2e2' + 3e|) - o. 

Let S be the solution set in M.^ to the above equation when = 1- By equation 
(39), it is clear that S is mirror symmetric about the ^1^2, CiCs and ^2^3 coordinate 
planes. Figure 1 below illustrates 5* in the quadrant ^1 > 0,^2 > 0,^3 > 0, and in 
this quadrant we see that S has one singular point ^sing G S. 




Surface S is defined implicitly by /(Cij '?2i "^s) = and singular points are charac- 
terised by V/ — 0. This yields ^sing — (y^iO; (For an alternative way to 
solve this point, see [Dah04, Lemma 4.2 (Hi)] ) Using computer algebra and the 
arguments used to prove Theorem 3.3 we may compute dim when Co = 1 and S 
intersects one of the coordinate planes {C,; = 0}f^]^. In these intersections we obtain 
dimV^ — 1 except at the singular point ^sing where dirnVj = 2. □ 
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4. Determining the medium from the Fresnel surface 

As described in the introduction, the new proof of imphcation (Hi) ^ (ii) in the next 
theorem is the first main result of this paper. Regarding the other imphcations let 
us make a few remarks. Implication (ii) (i) isa standard result for the Hodge star 
operator on a 4-manifold. The converse implication (i) (ii) is less well known. 
The result was first derived by Schonberg [Sch71]. For further derivations and 
discussions, see [Jad79, Rub02, HO03]. Below we will give yet another proof using 
computer algebra. The proof follows [HO03] and wc use a Schonbcrg-Urbantke-like 
formula (see equations (43)-(42)) to define a metric g from k. However, the below 
argument that g transforms as a tensor seems to be new. For a different argument, 
see [HO03, Section D.5.4]. 

When a general (2) -tensor k on a 4-manifold satisfies = —/Id as in condition (i) 
one says that n satisfies the closure condition. For physical motivation, sec [HO03, 
Section D.3.1]. Let us emphasise that Theorem 4.1 is a global result. The result 
gives criteria for the existence of a Lorentz metric on a 4-manifold. In general, we 
know that a connected manifold M has a Lorentz metric if and only if M is non- 
compact, or if M is compact and the Euler number xi^) is zero [MS08, Theorem 
2.4]. Let us also note that if J is an almost complex structure on a manifold A'/, 
that is, J is a (J)-tensor on A'/ with = — Id and dim A'/ > 2, then M is orientable 
[Hsi95, p. 77]. It does not seem to be known if the analogous result also holds for 
(2) -tensors, that is, if the closure condition on a 4-manifold implies orient ability. 

Theorem 4.1. Suppose N is an orientable A-manifold. If k & Vt2{N) satisfies 
= 0, then the following conditions are equivalent: 

(i) = —/Id for some function f G C°°{N) with f > 0. 
(ii) There exists a Lorentz metric g and a nonvanishing function f S C°°{N) 
such that 

(40) K = /*g. 

(Hi) ^^'k = and there exists a Lorentz metric g such that 

F(«) = 

where F{n) is the Fresnel surface for k and F{*g) is the Fresnel surface 
for the -tensor *g. 

Moreover, when equivalence holds, then metrics g in conditions (ii) and (Hi) are 
conformally related. 

Proof. For implication (i) => (ii) let rj = f~^/'^K whence jf' = — Id, and let h be 
an auxiliary positive definite Riemann metric on N . Let ^ be an atlas given by 
applying Lemma 4.2 to rj. For the local claim, let (C/, a;*) be a chart in £/' , and in 
this chart let 77 be represented by 3 x 3 matrices ^ and J(f . With computer algebra 
we then obtain 

(41) ^4,„(C,^,C,0 = sgndet^ (G»^^„6)', i^^\U), 
where G = {G"'^) is the 4x4 matrix 



(42) G = 



1 1 


/ det^ 




(det/i)V4 |det^|V2 




-^'J -F(det^)-^fc'fc^ ) 
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and /c* = s^^'' ^Sbcd-^'^'^ for i S {1, 2, 3}. Using a Shur complement [Pra94, Theorem 
3.1.1] we find that 

detG = - ^ 



dot h 



Hence det G < 0, so matrix G is invertible and has constant signature ( — h ++) or 
(H ) inU. Let Gij be the ijth entry of the inverse of G. In U we define 

(43) g = auGijdx'' ® dx\ 

where constant au S {—1, 1} is chosen such that g has signature ( — h ++). Then 
g defines a smooth symmetric (j) -tensor in U with signature ( — h ++), and by 
computer algebra we have 

(44) 'q\u = -sgndet.e/ *g . 

This completes the local claim in (i) => (ii). For the global claim, let (?7, .t') and 
(?7,x*) be overlapping charts in £^ , and in these charts let G"^^ and G^^ be defined 
as above. Since '^h,r] is a tensor, equation (41) implies that 

(45) sgndet^ (G^^"6^,)2 = sgndet ^ (ff^ ^ ^^rC.)' 

for all ^ = ^idx^ G A^(C/n{7). Since G"** is non-degenerate we can find a ^ such that 
the left hand side is non-zero. Thus sgn det ^ = sgn det j^/ inUCMJ and sgn det £/ 
in equation (44) defines a smooth function N ^ R. By Theorem 2.3 (in) (i) 
there exists a smooth nonvanishing function A : U (lU ^ R such that 

Equation (45) implies that function A can only take values {— 1,-|-1}. Thus 

Since auGij and cr^Gij both have signature ( — h -|— 1-). It follows that Xau = 
in [/ n [/, and equation (43) defines a tensor on A''. This completes the proof of 
implication (i) ^ (ii). 

Implication (ii) =J> (iii) follows by Propositions 2.2 and 3.6. 

For the proof of implication (iii) (i) we first establish two subclaims: 

Claim 1. The (Q)-tensor is pointwise proportional to ^^^^ by a non-zero 

constant. 

Let p G N. By Proposition 3.4 we only need to show that there exists a A e M such 
that 

%,*g{C,^,^,0 = A^^g,«(^,C,e,0, CeA^(Ar). 

Let a;* be coordinates around p such that g\p = A;diag(l, —1, —1, —1) for k G {±1}. 
In these coordinates, let J^gf^g and ^^g^^' be components for the symmetric (q)- 
tensors ^^g,*g|p and ^g^slp, so that 

%,*gii, = ^;%%^j^kii, e> a = ^fiii^ikii 

for ^ = ^jrfa:*|p. Using these components, let P, Q be the polynomials P, Q: — >■ M, 
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where S R, ^ = (6,6,6) € K^. By Proposition 3.6, 

= -(6-1^1)' (6 + 1^1)', 

for all (6,^) e when |^| is the Euclidean norm of ^. Thus P(l, 0,0,0) 7^ so 
dx'^lp i Fp{^g) = Fp{K) whence ^ 0. For each ^ e R^, (3(6, €) is then a 

fourth order polynomial in 6 with coefficients determined by ^ € M"^. Hence there 
exists continuous maps 

n-.R^^C, {1,2,3,4} 

so that for all ^ e M^, {ri{^)}j^-^ are the roots of Q{-,^) [NP94]. For each ^ e 
there exists a a(^) e K such that 

4 

(46) Q(6,0 = a(0 11(^0 -^*(^))' ^0^^- 

Applying d'^/d^Q to both sides implies that a(^) = k*"^- particular, the map 
4 1-^ a{^) is constant and non-zero. Let 11 = Since P and Q have the same 

zero set, there exists functions Sj : -> {—1, 1} such that 

n{^) = s,(OI€|, CeM^^e{l,2,3,4}. 

We know that ]R^\{0} is path connected. Hence ]R^\{0} is connected. For a con- 
tradiction, suppose that .s,;(M''\{0}) = {-1,+1} for some i e {1,2,3,4}. Then 
]R^\{0} = C/+ U U- for open, non-empty and disjoint sets U± defined as 

U± = {^ e M^{0} : ±r,(|) > 0}. 

It follows that there are constants si, S2, S3, S4 G {—1, +1} such that 

(47) r,(0 = s.lll, I eM^^e {1,2,3,4}. 

Let a be the number of Sj with Sj = 1. If ^ e M'^\{0}, then polynomial P{-,^) has 
two distinct roots ±|^|. Hence cr = or ct = 4 are not possible, so cr e {1, 2, 3} and 
by equation (46), 

Q(6,0 = M(6-Ill)"(6 + I€l)'"" 

f/i (6^ + 26141(6' -l^n), if^ = i, 

(6 -141^- 26141(6' -l^n), if^ = 3, 

for all (6,4) £ H^^- Since Q is a polynomial, we know that t 1— ?• (3(1, 0,0) is 
smooth near 0. This is only possible when a = 2, and Claim 1 follows. 

Claim 2. At each p G N there exists a non-zero A e M such that k\p = X*g \p. 

Let p G N. By Proposition 2.4 (in) there are coordinates around p such that 
g\p = A;diag(l, —1, —1, —1) for some k € {±1} and s^\p is diagonal. For 4 = ^idx^\p 
we then have 

(6 6 6 6 = ^;i%CjS.k^i, %At 6 6 6 = ^^;i''6666, 

where and ^'g^^^^ are components for ^g,*^ and ^g^^. in coordinates a;*. By 

Claim 1 there exists a A e R\{0} such that 

%fg = X^;f, 0<i<j<k<l<3. 

Moreover, ^^^k = 0. We then have 36 polynomial equations for k. Using the Grobner 
basis (see Appendix A) for these equations we find that the equations have a unique 



16 



DAHL 



real solution for n and this solution is given by k\p = A *g \p. This completes 
the proof of Claim 2. 

By Claim 2, there exists a map X: N such that k = X*g whence = 

— A^Id. To see that is smooth it suffices to note that A^ = —g trace k^. This 

completes the proof of implication (in) (i). 

When equivalence holds, Proposition 3.6 and Theorem 2.3 imply that the Lorentz 
metrics in conditions (ii) and (Hi) are conformally related. □ 

The next lemma was used to prove implication (i) (ii) in Theorem 4.1. In the 
proof of the lemma, Claim 1 is based on [HO03, Sections D.4-D.5]. 

Lemma 4.2. Suppose N is an orientable i-manifold and k G 02(-^)- V ^o,s 

no skewon component and = — Id, then N has an oriented atlas .'7 with the 
following property: Each p S N can be covered with a connected chart (C/, a;') € £^ 
such that if £/, ^, ^ represent k in U, then 

(i) is invertible in U. 

(ii) In U there exists a smoothly varying antisymmetric 3x3 matrix ^ such 
that 



Proof. Let us first make an observations: Suppose {a;'}f_Q are arbitrary coordinates 
for N and £/, SS, ^€ , '2 arc 3x3 matrices that represent k in these coordinates. 
Then Proposition 2.4 (i) implies that k'^ = — Id is equivalent with 

(48) '£'^+s^S§ = -Id, 

(49) W + = 0, 

(50) '^j^ + jaf'^'^ = 0. 

Let ^ is a maximal oriented atlas for A''. The proof is divided into two subclaims, 
Claim 1 and Claim 2. 

Claim 1. For each p G N there exists a connected chart (t/, x*) that satisfy 
condition (i) and there exists a chart {W,y^) G .% with U (IW ^ such that the 
transition map n- is orientation preserving. 

By Proposition 2.4 (ii) we can find a connected chart {U,x^) that contains p and 
where matrix £/ for k is diagonal at p. The rest of Claim 1 is divided into four 
cases depending on the eigenvalues of s^\p. 

Case A. Suppose all three eigenvalues of £/\p are non-zero. Since eigenvalues 
depend continuously on the matrix entries [NP94], we can shrink U and part (i) 
follows. Claim 1 follows by possibly reflecting the a;^-coordinate. 

Case B. Suppose .s/\p has two non-zero eigenvalues. By permutating the coor- 
dinates (see equation (21)) we may assume that s^\p = diag(ai, a2, 0) for some 
0-1,0,2 7^ 0. Writing out equations (48)-(50) with computer algebra gives 

75 1 — ^ 2 — ^ 1— ^ 2 — U, 3) — — i 

at p. The last equation contradicts that ^ is real. Case B is therefore not possible. 

Case C. Suppose s^\p has one non-zero eigenvalue. As in Case B, we can find a 
chart {U,x^) for which s^\p = diag(ai,0,0) for some ai ^ 0. Writing out equations 



DETERMINATION OF ELECTROMAGNETIC MEDIUM FROM THE FRESNEL SURFACE 17 



(48)-(50) as in Case B gives 

^1 1 = ^ = <^3 1 = 0, ^11 9^ 0, 7^ 0, 

at p. Let {x^}i^o be coordinates around p defined as 

3j — X I X ^ 

F = x\ ie{l,2,3}. 

In these coordinates, matrix £/\p has determinant —=^11(^^2)^) which is non-zero, 
and Claim 1 follows as in Case A. 

Case D. Suppose all eigenvalues of £/\p are zero. Then £/\p — and equation (48) 
implies that (det'^lj,)^ = —1. This contradicts that "^jp is a real matrix. Case D 
is therefore not possible. 

Claim 2. Lc!t be the collection of all charts (f/, a;*) as in Step 1 when p ranges 
over all points in N. Then satisfies the sought properties. 

Let {U,x'^) and (?7,.x*) be overlapping charts in 3^. By Claim 1 and [Tro94, Lemma 
13.9], each chart U and U is compatible with all charts in J^q. Hence the transition 
map a;* x^ is orientation preserving, and ^ is oriented. We know that each 
chart in ^ satisfies property (i), and property (ii) follows by defining ^ = "^jz/. 
Indeed, J^^ is antisymmetric by equation (50), and the expression for follows by 
equation (48). □ 

5. NON-INJECTIVITY RESULTS 

Implication (Hi) (ii) in Theorem 4.1 shows that for a special class of medium, 
the Fresnel surface determines the medium up to a conformal factor. In this section 
we will describe results and examples where the opposite is true. In the below we 
will see that there are various non-uniquenesses that prevents us from determining 
K (or even the conformal class of k) from only the Fresnel surface F{k). 

Let us study the non-injectivity of the two maps in the diagram below; 

(51) K ^(k) F{k), k e nl{N), 
where ^^(k;) is the Tamm-Rubilar (Q)-tensor density induced by k. 

5.1. Non-injectivity of leftmost map. Let us first study the non-injectivity of 

the leftmost map in diagram (51), that is, the map 

(52) K <^{k), Kenl{N). 

Parts (ii) (iv ) in the next theorem describe three invariances that make the map in 
(52) non-injective. The first two parts are well known [HO03, Section 2.2]. However, 
let us make three remarks regarding part (iv). First, an interpretation of part (iv) 
is as follows: If F, G solve the sourceless Maxwell equations in medium k, then G, F 
solve the sourceless Maxwell equations in medium . In this setting, part (iv) 
states that both media have the same Fresnel surfaces. Second, suppose *g is the 
(2)-tensor induced by a pseudo-Riemann metric g. Then *2 = ±ld, so = ±*g, 
whence '^{*g) and $f(*~^) are always conformally related. Part (iv) states that this 
is not only a result for Hodge star operators, but a general result for all (2) -tensors. 
Third, the proof of part (iv) is based on computer algebra. Of all the proofs in this 
paper, this computation is algebraicly most involved. For example, if we write out 
equation (53) as a text string, it requires almost 13 megabytes of memory. 

Theorem 5.1. Suppose k G fl'2{N) where N is a A-manifold. Then 
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(i) c^IJk) = /^^(k) for all f G C°^{N), 

(ii) = 0, 

(ni) ^(k) = (k + / Id) for all f G C^iN), 

(iv) ^{k^^) = ^#'(— (det «;)~^/^k) when k is invertible. 



Proof. Part (i) follows by the definition, and parts (ii) (Hi) are proven in [HO03, 
Section 2.2]. Therefore we only need to prove part (iv). Let adj k = detKK"^ be 
the adjugate of n. By part (i) \t suffices to show that 

(53) (dct k)2^^^'=' + = 0, 0<z<j </c</<3, 

where and are components of the Tamm-Rubilar tensor densities of n 

and adj k, respectively. The motivation for rewriting the claim as in equation (53) is 
that now both terms are polynomials. By using the method described in Appendix 
B we obtain that equations (53) hold, and part (iv) follows. □ 

Theorem 5.1 (ii) shows that if we restrict the map in equation (52) to purely skewon 

tensors, wc do not obtain an injection. The next example shows that the same map 
is neither an injection when restricted to tensors of purely principal type. 

Example 5.2. On iV = M x with coordinates {x*}?^o> ^ be the (2)-tensor 
defined by 3 x 3-matrices 



^ = 03x3, ^=10 AaL'T^I -2-1/3 As I , ^ = 

where parameters Ai, . . . , A5 € R are arbitrary. Then n has only a principal part, 
det/t = 1, and 

for any pseudo-Riemann metric h on N. □ 

When proving implication (Hi) => (i) in Theorem 4.1 we need to assume that n 
has real coefficients. In fact, for (2) -tensors with complex coefficients a decompo- 
sition into principal-, skewon-, and axion components does not seem to have been 
developed. The next example shows that there are non-trivial complex tensors 
whose Fresnel surface everywhere coincides with the Fresnel surface for the stan- 
dard Minkowski metric go = diag(l, —1, —1, —1). (For k G OjC-^j^) define the 
Fresnel surface using the same formulas as for real k.) 

Example 5.3. On = M x M'^ with coordinates {x^}f^Q, let k be the (2)-tensor 
with complex coefficients defined by 3 x 3-matrices 
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where z is an arbitrary function z: N ^ C\{0} and i is the complex unit. At each 
p G N the Fresnel surface is then determined by 
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where £,idx^ S Ap{N), and 

(1 + 6z^f (5 - 126z3 + 684^6 _ 6482^) 

trace «; = 0, det k = — ^ r^, . 

46656 

Prom the latter equation we see that for specific values of z, tensor k can be non- 
invertible as a linear map □ 

5.2. Non-injectivity of rightmost map. The next example shows that there are 

skewon-frcc (2) "tensors ki and K2 that have the same Fresnel surfaces, but their 
Tamm-Rubilar tensors are not proportional to each other. This shows that the 
rightmost map in equation (51) is not injective. Let us point out that this contra- 
dicts the first proposition in [PSW09] whose proof does not analyse multiplicities 

of roots to the Frcsncl equation. 

Example 5.4. On iV = M x with coordinates {x'jf^Q, let ki be the (2) -tensor 
defined by 3 x 3-matrices 



2^ 1 1 , ^1 = < 




2 





For the Euclidean metric go on N wc then have 

%o.nAi,ti,o = (eo-ei)(co-e2)', CeA^iV). 

To exchange the role of ^1 and ^2, we perform a coordinate change xq i— )■ xq, 
xi i-> X2, X2^ xi, 2:3 ^ X3. With this as motivation we define K2 as the (2)-tensor 
defined by 3 x 3-matriccs 

0' I , = I I , ^2 = '^2^ 



2 


1 













-1 


1 




Then 



Here ki and K2 arc not proportional, their Tamm-Rubilar tensor densities are not 

proportional, but their Frcsncl surfaces coincide. 

Both Ki and K2 have 1 has an eigenvalue of algebraic multiplicity 6. Hence 

det Ki = det K2 = 1, trace ki — trace K2 = 6, 
and for the trace-free components Ki = Ki — Id we have det Kj = 0. □ 
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Appendix A. Grobner bases 

To prove implication (in) (i) in Theorem 4.1 wc use a Grobner basis to solve 
a system of polynomial equations. In this appendix we collect the results about 
Grobner bases that are needed for this step in the proof. These results are standard 
and can, for example, be found in [CLO07, pp. 5, 29-32, 76-77]. 
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By C[a;] we denote the ring of polynomials with complex coefficient that depend on 
variables xi,. . . , xn G C where A'' > 1. A (polynomial) ideal is a subset / C C[x] 
such that 

(t) e /, 

(a) f,g & I implies that f + g & I, 
(Hi) f & I and g G C[x] implies that fgGl. 

Theorem A.l (Hilbert basis theorem). If I C C[x] is an ideal, then there exists 
finitely many polynomials /i, • • • , /« € / such that 

(54) I = eC[x]|. 

When Theorem A.l holds, we say that polynomials f\,. . . ,fs form a basis for ideal 
/. Conversely, if /i, . . . , fs arc any polynomials in C[x], the ideal on the right hand 
side in equation (54) is denoted by (/i, . . . , fg) and called the ideal generated by 
polynomials /i, . . . , /«. The affine variety defined by polynomials /i, . . . , e C[a;] 
is the subset V{fi,. . . , /^) C C^, 

V{h,...,fs) = {xgC" :Mx) = --- = Mx)=0}. 

The next proposition gives sufficient conditions for two systems of polynomial equa- 
tions to have the same solution sets. 

Proposition A. 2. Suppose /i, . . . , and gi,. . . ,gt are polynomials in C[x]. If 

{fi} and {gi} generate the same ideal, that is, 

(/i,---,/s) = {gi,---,gt), 

then 

V{fi,...,fs) = V{g,,...,gt). 

We will not give a precise definition for a Grobner basis. The key property for 

Grobncr bases is collected in the next proposition. 

Proposition A. 3. Let /i,...,/^ G C[x] be polynomials such that (/i, •••,/«) 7^ 
{0}. Then there exists polynomials gi,. . . ,gt G C[x] such that 

(/i,---,/s) = {gi,...,gt)- 
Polynomials gi are called a Grobner basis for the ideal (/i, . . . , 

Even if computation of Grobner basis computation is supported by modern com- 
puter algebra systems, their computation can in practice be very time consuming. 
The motivation for using a Grobner bases is that they typically simplify the solu- 
tion process for polynomial equations. Thus one can think of Grobner bases as a 
way to simplify polynomial equations without changing their solution set. This is 
illustrated in the next example. 

Example A.4. Let 5 C be all (x, y, z) e such that 

(55) xyz = 1, xz'^ = y^, = xy. 

By elementary manipulation, we see that S = {(1,1,1)}. To illustrate how to 
determine S using a Grobner a basis, let us first note that S = V{fi,f2, fs) n M^, 
where 

fi=xyz-l, f2=xz'^-y'^, fs = z^-xy. 
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With computer algebra we find that a Grobner basis for (/i, /2, /s) is given by 

Propositions A. 2 and A. 3 imply that V{fi,f2, /a) = V{gi,g2,g3)- Hence 5" coincide 
with real solutions to polynomial equations 

(56) z^ = l, y^ = z, x = y''z. 

These last equations are easily solved and we find that S = {(1,1,1)}. 

If we compare the original equations (55) to equations (56) computed by a Grobner 
basis, we see that the latter ones are more easier to solve since they can be solved 
by backsubstitution. □ 



Appendix B. Verifying very large polynomial identities 



The proof of Theorem 5.1 (iv) reduces to proving equations (53) which consists of 
35 polynomial identities in 36 variables. If we write out these polynomial identities 
as text strings, they occupy almost 13 megabytes of memory. Due to this size, 
Mathematica (version 7.0.1) was not able to verify the identities in a reasonable 
time. In this appendix we describe a recursive method that is able to verify these 
identities. On a computer with two Intel E8400 3GHz processors and 3.7 gigabytes 
of RAM the method finished in 10 hours. The method relies on the following 
corollary to Taylor's theorem with a Lagrange error term. 

Proposition B.l. Suppose f is a polynomial f : -^R in variables xi,. . . , xn € 
R. Furthermore, suppose that 

(i) There exists a finite G {1, 2, . . .} such that 



■j^{x) = for all X e 



nJV 



dx^ 

(ii) Polynomials Z°,..., Z^'^ : R^'^ R defined as 

Z'{y) = f^(0,y), 2/eK^-\ r = Q,...,K-l 

are zero polynomials. 
Then f is the zero polynomial. 

Proposition B.l shows that to verify identity / = we only need to verify identities 
Z° = 0, . . ., Z^~^ = 0. Since these identities are obtained by differentiating / and 
by setting one variable to zero, they are typically shorter and easier to manipulate 
than the original /. By recursively applying Proposition B.l, the proof of / = 
divides into smaller and smaller polynomial identities that eventually can be verified 
using Mathematica's internal Expand routine. The implementation details are as 
follows. To verify f — we applied Proposition B.l recursively until the polynomial 
had less than 27 variables (out of 36 original) . For each application of Proposition 
B.l we used the (non-optimal) constant K = 5. 
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